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Abstract 

A comprehensive number of one-loop integrals in a theory with Wilson fermions at 
r = 1 is computed using the Burgio-Caracciolo-Pelissetto algorithm. With the use of 
these results, the fermionic propagator in the coordinate representation is evaluated, thus 
making it possible to extend the Liischer-Weisz procedure for two-loop integrals to the 
fermionic case. Computations are performed with FORM and REDUCE packages. 
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1 Introduction 



Perturbative calculations in lattice gauge theories (for a review, see [T]) are of interest from 
several points of view. 

Firstly, they are needed to determine the Alat parameter of QCD in the lattice regularization 
and its relation to the respective value Aqcd in the continuum theory 

Secondly, every lattice action defines a specific regularization scheme, and thus one needs a 
complete set of renormalization computations in order for the results obtained in Monte Carlo 
simulations be understood properly. Perturbation theory is required to establish the connection 
of the matrix elements computed on a lattice with their values in the continuum theory [2J, [3]. 
In this connection, it should be emphasized that the use of one-loop perturbative renormalization 
constants gives rise to large systematic uncertainties in lattice calculations of the momenta of 
hadronic structure functions [3J and respective two-loop computations are needed. 

Thirdly, perturbative calculations provide the only possibility for an analytical control over 
the continuum limit in QCD. One can also mention anomalies, proof of renormalizability, 
Symanzik improvement program and other fields of application of lattice perturbation theory. 

Here we consider one- and two- loop diagrams with Wilson (r = 1) fermions at zero external 
momenta jl]. We outline the Burgio-Caracciolo-Pelissetto (BCP) method j5] of calculations of 
one- loop integrals and describe the respective computer algorithm [B]. This algorithm allows to 
compute the fermionic propagator in the coordinate representation and, therefore, to extend the 
Liischer-Weisz (LW) method [7] to the fermionic case; such extension is presented in Section 4. 



1.1 Notation 

We use the following designations: h stands for the set ni, n 2 , n 3 , n 4 ; x = (x\,X2, £3, £4), where 
x^ are integer- valued coordinates of an infinite four- dimensional lattice A = {x : x^ G we 
also need the lattice A' = A\{0} with removed site x=(0, 0,0,0); 

\x\ = \x X \ + \x 2 \ + \x 3 \ + \x 4 \, [x n ] = X\ + X% + £3 + x\ . (1) 

Then we give the expressions for the denominators of bosonic and fermionic propagators, 

As{k) = 4 + fj? R — cos(/ci) — cos(A;2) — cos(k 3 ) — cos(k i ); (2) 
4 

A F (k) = 10 - 4 cos (M + cos(fe /i ) cos(fc„) + fi 2 R 

where \l r is the fictitious mass for infrared regularization. We also use Dp = 2 Ap and 
Db = 2 Ab normalized in the standard way {Db{f){K) — 1/k 2 as k ^ 00). These propagators 
in the coordinate representation are defined as follows: 

G B (F)(x) = / TTT^T^ 77T, (3) 

JBZ (27T 4 ) D B {F)[k) 



where BZ is the Brillouin zone, BZ = <p : < £>u < — fj 

1 a a> 



p M = -sinf^V p 2 = Y Pi; (4) 



mu=l 



a \ q 

a is the lattice size. In this work we set a = 1 for the sake of simplicity. 
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2 The Burgio— Caracciolo— Pelissetto method 
2.1 Bosonic Intefgrals 

The integrals under study are denned as follows: F(q; h) = lim^o Fs(q; n), where 

jp / ~ % f „ cos(A; 1 )" 1 cos(A; 2 )" 2 cos(A;3) na cos(A;4) ri4 
^ n) = 7bz A^ 3 ' ® 

Here 5 is an infinitesimalparameter for an intermediate regularization [5]. This parameter 
makes it possible to derivqj the recursion relations of the form 

F{q;...,n^...) = F{q; ...,n^ - 2, ...) - (6) 
(n M - l)F(g-l;...,n M -l,...) | (n^ - 2)F(g - 1; ...,n^ - 3, ...) ^ > ^ 



g - 1 + 5 q-l + S 

With these relations and similar relations for n M < 1, one can express the integrals ([51) in terms 
of the quantities 

G5(g '^ )= L(^(A^' (7) 

Up to terms of the order 0(n 2 R ) and 0(5), this expression has the form 

Fs(q,n)= £ 4; ) (^^)^(r,0) + E4r ) (^^)^(r,^), (8) 

r=g— ni— ri2— ri3— n4 r=l 

where A^(5,h) have a pole singularity in 5, and are polynomials in n 2 R . As for the 

function G$(r,fi R ), the domains r > and r < should be considered separately. At r > 0, 6 
can be safely set to zero and the function Gs(r,fj, R ) should be expanded in powers of fi R 2 : 



Gs{r, /4) 



1 



r-2 



+ J(r) + 0(/i|) + O(<5), (9) 



(2vr) 2 r(r) 

where b n are the coefficients of the asymptotic expansion at z — > oo of the functional 

exp(-4,)/ * W . + 1 + l + (10) 

Zc* = hi(yU^) +C, and C is the Euler-Mascheroni constant. At r < 0, /jlr can be safely set to zero 
and the function G$(r, 0) should be expanded in 5 as follows: Gg(r, 0) = B{r) + J(r)5 + 0{5 2 ). 
The functions J(q), in their turn, obey recursion relations of the type 

c (q)J(q) + ci(g) J{q + 1) + c 2 (q)J(q + 2) + c 3 (q)J(q + 3) + c 4 (g) J(q + 4) = (11) 

derived in [5J; the explicit expressions for the coefficients Cj(g) can be found in [B]. Thus we 
express J(g) at g > 4 and at q < in terms of J(0), J (I), J (2) and J(3). It should be noted 
that J(0) does not appear in ultimate expressions for the integrals (J5j. Then one can introduce 
the values 

Z ='^, F = 4vr 2 J(2), Z x = 32 J(3) - 8J(2) + + - (12) 
which are equal to P Z « 0.15493339023, Z x « 0.10778131354, F « 4.369225233874758. 



1 Using integration by parts 
2 Io(z) is the Infeld function. 
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2.2 Fermion Integrals 

In the fermionic case, we consider the quantities F(p, q; h) = lim^o F$(p, Q] n), where 

. , r d 4 k cos ni (k 1 )cos n2 (k 2 )cos n3 (k 3 )cos ni (k 4 ) . . 

Fs(p,q;n) = hm / - — — . (13) 

With the recursion relations similar to (JBJ), these integrals are expressed in terms of the functions 
r d 4 k 1 

Gs(p, q) = — tj, which can be represented in the form 

J (27TJ /\ B /\p 

G s (p,q) = D(p,q;^ R ) + B(p,q) + 5(L(p,q;fi 2 R ) + J(p,q)) + 0(5 2 ), p < 0; (14) 
G s (p,q) = D(p,q-fi 2 R ) + J(p,q) + 0(6), p>0. 

The divergent parts D(p,q;fM R ) and L(p,q;fi 2 R ) in the domain of interest can be determined 
by a straightforward procedure [5], whereas the functions B(p,q) and J(p,q) obey recursion 
relations of several types. These relations and the procedure of their derivation were described 
in [5] ; their explicit form (very cumbersome) is given in [B] . With the use of these relations, the 
functions F(p,q;n) can be represented (see [T], [S]) as linear combinations of the constants F , 
Z , Zi and 

y ° = ^-i^ ^sri*-^- 1 -^* - 1 ^^ (15) 

Y, = i-i-z„-ij(-l,2) + ij(0,l)-i 7(1,-2) 

17 1 25 

__ J(i, 0) + 4 J(l, 1) - - J(2, -2) + - J(2, -1) - 4 J(2, 0), 

Y 3 = l — -F + — Z - — J(-l,3) + — J(0,1) + ^- J(0,2) + — J(l,l); 

384vr 2 128tt 2 96 48 V ' ; 192 v ' ; 48 V ' ; 48 V ' ; ' 

n = ^y^, F 5 = J(1,-1), F 6 = 2J(l,-2), Y 7 = ^^l, (16) 

F 8 = J(2,-2), Y 9 = J(3 2 ~ 2) , Yi = J(3,-3), r n = 2J(3,-4). 

The respective codes can be found on the web page of the ITEP Lattice group 
http://www.lattice.itep.ru/~pbaivid/lattpt/, The results stored there are as follows: 

(i) the program for a computation of F(p, q; h) at < p, q < 9 and n\ + n 2 + n 3 + n 4 < 25; 

(ii) the values of the functions J(p, q) and B(p, q) at —26 < p < 0, —56 — 2p < q < 34 and 
the values of J(p, q) at 1 < p < 9, —28 < q < 33 — p; and (iii) The explicit expressions for 
F(p, q; h) at some particular values of p and q and all n\ < 6. 

3 The Liischer— Weisz method 

To outline the LW method [7] of computation of two-loop diagrams in the coordinate represen- 
tation, we consider the diagram in Fig.l, given by the expression 
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B 



X . 



(17) 



X 




In the bosonic case, Liischer and Weisz calculated 
Ab(0) and its asymptotic expansion when p — > 0; 
they used the following representation: 



A B (0) = GU0)+J2G 3 a 



[x 



xeA' 

+ E (g|(x)-gL(x) 
+ E - 

where J-jv = {x : + |x 2 | + |a?3 1 + \xi\ < N}, and Gr as (:r) is an asymptotic approximation of 
Gb{x) when x — >■ oo, 



Figure 1 



+ 



'2[x 4 ] - [x : 



212" 



„412 



214 



+ 40 



„217 



+ 16 



215 



48 



,216 



213 



(19) 



In the domain J 7 ^, the propagator Gb{x) can computed by the recursion formulas 



G B {x + ft) = G B {x - ft) 



2x, 



^2u=l X v ) A=l 



E(G , B(x)-G B (a;-A)), 



(20) 



which allow to express it in terms of Gb(0, 0, 0, 0) = Z and G B (1, 1, 0, 0) = — 1/4 + Z\ + Z . 
The domain {Ayj 7 ^} is chosen so that the propagator is fitted by its asymptotic expression 
( IT9l) with a sufficient precision making it possible to neglect the third sum in the formula ( fl"8l) . 
Then the first sum can be calculated exactly using the summation formulas derived in [7] and 
the second sum can be expressed in terms of Z\ and Z 2 by employing the relations ( 120]) . 
It should be emphasized that Ab{0) is the coefficient of the expansion 

A B {p) = \{A B (0) + (pa) 2 U + B 1 hx(pa) 2 } + 0((pa) 4 )), 

where a is the length of a link; that is, A B (0) is the coefficient of the divergent part. It 
does not vanish, though the "bigmac" diagram depicted in Fig.l converges in the dimensional 
regularization provided that all masses and the external momentum are equal to zero. 



4 Two-loop fermionic integrals. 

In the fermionic case, calculations are performed by the same procedure, however, we have no 
recursion relations similar to ( |20j) . The fermionic propagator in x- representation 

_ , . f d A k cos(fciXi) cos(fc 2 a;2) cos(fc 3 x 3 ) cos(fc 4 X4) . . 

G F ( Xl , x 2 , x s , x 4 ) = J — — (21) 

is expressed in terms of the quantities 

f d 4 k cos" 1 (A; 1 )cos n2 (A;2)cos ri3 (fc3)cos n4 (A;4) , , 

F(p, q; m, m, n 3 , n 4 ) = J A^A^ ( ^ 



by making use of the relations 

[n/2]— 1 / -, \k+l 

cos(nx)=2"- 1 cos"x + - £ ^f— Ct fc _ 2 (2cosxr 2fc ~ 2 . (23) 

To employ the LW method outlined above, we compile a table of values of Gp{x) over the 
domain X\ > x 2 > x 3 > x 4 > 0, |x| < 48 and derive an asymptotic approximation of Gp(x) at 
\x\ — > oo up to the terms of the order l/[x 2 ] 4 . To treat integrals with nontrivial numerators, we 
should also compile the tables of the values 




f dp (e~ ipx — 1) 

Kb ^ = 1W)~b^W LBm = }W)~ D% [F]{P { ' (24) 

Each of these tables involves 14147 entries, each entry is a linear combination of the constants 

F , Z , Zi, Y , Yi, ...lii, t — \n > an d 1 with rational coefficients; from 5 to 20 MB per table in 

(2tt) 2 

size. Fortunately, they can be conveniently treated with FORM [8]. 

The precision of 20 significant digits in determination of the constants Y± Y\\ [T], jS] is 
not sufficient for computation of Gp(x) at \x\ > 6. Using the procedure proposed in [7] for 
calculation of Z and Zi, we obtain 

Y A = 0.08539036359532067913516702888533412058194147127443265(1 

Y 5 = 0.46936331002699614475347539705751803482046295887523184(1 

Y 6 = 3.39456907367713000586008689702374496453685272313733503(1 

Y 7 = 0.05188019503901136636490228766471579940968012757291508(1 

Y 8 = 0.23874773756341478520233613930386970445280194983477988(1 

Y 9 = 0.03447644143803223145396188144243193600121277124715784(1 

Y 10 = 0.13202727122781293085314731098196596971197144795959477(1 

Fn = 0.75167199030295682253543148590778110991011277193144803(1 

At \x\ > 48, Gf(x) is approximated by the function 

<*■»« - p] + + H ff - » ft + » ft _ ^) + (25) 

To provide an example, let us consider the following two-loop fermionic integrals: 

Q BBB = f d*k A 4 klql 

Vl Jbz (2tt) 4 (2tt) 4 ^ D s{k)D b{q)D b{x) 1 J 



(2vr) 4 (2vr) 4 ^ D B (k)D B (q)D F (r) 



and similar quantities with other combinations of bosonic and fermionic propagators. We can 
also consider 

d'k d\ * kie/i 



qBBF _ / ^ M (27) 

^ 2 7bz (27T) 4 (2vr)4 ^ D B (k)D B (q)D F (r) 1 J 
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etc. The results of the computations are as follows: 



Qf BB = 



-- 0.042306368(1) 
= 0.024555253(3) 
= 0.01173224(1) 
0.05462397818(1) 
0.0130373237(1) 



oi 

Q{ 
Of 
Ql 
Ql 



,FBB 



0.020079702(3) 

0.00969896(1) 

0.00576013(3) 

0.02659175158(3) 

0.0064945681(3) 



(28) 




FFB 



,FFF 



BBF 



FFF 



5 Summary and Outlook 

The BCP algorithm has been realized on a computer. The basic fermionic integrals G(p, q) are 
found over a sufficiently large domain of values of p and q. This allows (i) to express one-loop 
intergals involving fermionic denominators in terms of the constants F , Z , Z\ and Yq-Z-Yh and 
(ii) to compute Gp(x) at |x| < 96. 

The LW method is extended to the case of fermions; asymptotic behavior of the fermionic 
propagator at \x\ — > oo is found. Therewith, Gp(x) is expressed at \x\ < 48 in terms of the 
constants Y^-^Yu, the values of which are computed to a precision of 54 significant digits. This 
is really needed for calculation of two-loop integrals. A new feature of FORM - a possibility 
to work with database-like structures - proved to be useful for summation over the domain 
\x\ < 48. As an illustration, several two- loop fermionic integrals are evaluated at zero external 
momentum. 

Operations with a table of precise values of the functions Gb(f)(x) K b ^{x) and L B ^(x) al- 
low to compute one-loop and two-loop diagrams of the propagator type at nonvanishing external 
momentum. The work is in progress! 
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